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K ^ 7T7T Decay Amplitude on the Lattice 

N. Ishizuka ^ 

^ Center for Computational Physics, University of Tsukuba, Tsukuba, Ibaraki 305-8577, Japan 

Recent theoretical and numerical progresses of the lattice calculations of if — > tttt decay amplitude are reviewed. 



1. Introduction 

Despite the full understanding of the funda- 
mental theory of weak interactions, the non- 
leptonic decay of hadrons still remains as the least 
understood of weak processes, the most notable 
problems being the A/ = 1/2 rule and the value 
of e'/e. The predicament originates from the dif- 
ficulty of evaluating the K — > tttt decay ampli- 
tudes. In spite of recent progress of computa- 
tional power, the calculation of the decay ampli- 
tude by lattice QCD simulations is much more 
difficult than those of hadron masses. There are 
two difficulties in the calculations, one of which 
is the operator matching of 4-fermion weak oper- 
ators, and the other is the problems first pointed 
out by Maiani and Testa 

For the A/ =1/2 decay process, operator mix- 
ings with lower dimensional operators can occur. 
Since these operators do not contribute to the 
physical decay amplitude, we have to subtract 
the effect of these in a non-perturbative manner. 
The non-perturbative subtraction is not easy nu- 
merically, but explicit methods have been pro- 
posed and carried out in actual lattice calcula- 
tions. The details are reviewed in Ref. j^. A re- 
cent progress is the realization that twisted mass 
QCD approach may provide a powerful tool for 
solving the operator mixing 1^,^. 

The most serious problem in the evaluation of 
the decay amplitude was pointed out by Maiani 
and Testa |]l|. We have to calculate the de- 
cay amplitude into a two-pion state with non- 
zero momenta. However, it is non-trivial to ex- 
tract it from the corresponding Green's function 
K{0) —>■ 7r(p)7r(— p), since the two-pion state with 
non-zero momenta is an excited state in the two- 
pion system. This problem of extraction of the 



decay amplitude is one of the problems pointed 
out by Maini and Testa, which we shall refer to 
as "MT-1" in this article. 

To avoid this problem, in most of lattice cal- 
culations, the decay amplitudes are calculated 
either at an un-physical kinematics or through 
the K n amplitudes and chiral perturbation 
theory (CHPT) relating these amplitudes to the 
physical decay amplitudes. Using such an effec- 
tive theory, however, is the cause for large un- 
certainties of the lattice prediction of the decay 
amplitude. 

Even if we succeed in overcoming the MT-1 and 
extract the decay amplitude at the physical kine- 
matics from the Green's functions, another prob- 
lem exists. Since lattice simulations are carried 
out on a finite Euclidean space-time, the decay 
amplitude on the lattice does not directly give 
the physical one in the infinite volume Minkowski 
space-time. This problem of the relation between 
the lattice and the physical amplitude is another 
part of the Maini-Testa problem, which we shall 
call as "MT-2" in this article. 

Recently Lellouch and Liischer derived a re- 
lation between the two amplitudes Q. Their 
derivation does not rely on effective theories. Lin 
et.al. reached the same relation from a different 
approach, and extended it to general kinemat- 
ics 1^ . They also examined the conditions for the 
validity of the relation. 

In this article I focus on the Maini-Testa prob- 
lems. Recent theoretical and numerical progress 
for avoiding or solving the MT-1 will be discussed. 
We also show a brief derivation of the relation by 
Lellouch and Liischer. We refer to the reviews of 
weak matrix elements in recent lattice conferences 
by Lellouch § and Martinelli [|, and Ref. [| for 
other recent theoretical and numerical progresses. 
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2. The MT-1 

The MT-1 is the difficulty of extraction of the 
decay amplitude of the two-pion state with non- 
zero momenta. Originally Miani and Testa found 
the problem for the K — > tttt Green's function in 
infinite volume Euclidean space-time 0| . Here we 
consider a finite- volume Green's function given by 

Gn{t) = (0| (w)„(t) 0(0) K{tK) |0) , (1) 

where (7r7r)„(t) is the interpolating field for the 
n-th excited two-pion state given by (7r7r)„(i) = 
7r(p„,t)7r(-p„,t) and pi = {^Ti/Lf ■ n, 0(0) is 
the weak operator and K{tK) is K meson filed. 

Inserting the energy eigenstates, the Green's 
function (^ can be rewritten as 

Gn (t) = ■ • ■ (^^) ' (2) 

for the time region f 3> 3> i_R-, where 
= (i^j\0\K), Vnj = (0|(7r7r)„|(^^.), and 

Aj{t) — exp{—Ej ■ t). The state |(7r7r)^.) is the 
j-th two-pion energy eigenstate with the energy 
Ej, which is different from twice the pion en- 
ergy Ej due to the two-pion interaction on the 
finite volume. The energy Ej satisfies the Liischer 
quantization condition The state \K) is the 
zero momentum K meson state and G^^ {t) — 

{K\K\Q)-eM-mK-tK)- 

The crucial point is that Vnj 9^ 6nj in (||) gener- 
ally, because the interpolating field (7r7r)„ has no 
definite energy and it can emit the state |(7r7r)j) 
with j ^ n. Thus, Gn{t) contains many exponen- 
tial terms generally. The extraction of the decay 
amplitude from such a multi-exponential function 
is not trivial. For n = 2, for example, the dom- 
inant contributions of Gn{t) for the large time 
region come from the states with j = 0, 1, 2, and 
the effects from these states disturb extraction of 
the amplitude A special case is n = 0, for 
which we can extract the amplitude by a single 
exponential fitting, because contaminations from 
lower states are absent. 

3. Method for avoiding the MT-1 

There are three methods known for avoid- 
ing the MT-1. These are to calculate the 
Green's functions for the amplitudes of (1) K 



7r(0)7r(0) 0, (2) if TT 0, and (3) ^ 
tt(P)iM) at p2 ^ (27r/L)2 (SPQR kinemat- 
ics) |1^. As we noted in Sec. g, the MT-1 does 
not apply to the zero momentum two-pion. Fur- 
ther this is also true in case (3), because only 
one state is dominant for the large time region. 
Results of lattice calculations by the above three 
methods are reviewed in Sec ^. 

In the methods (1) and (2), chiral perturbation 
theory (CHPT) is used to reconstruct the physical 
decay amplitudes from those on the lattice. Here 
we consider, as an example, the case of operators 
for (27l,1_r) and (8l,1_r) representations after 
subtraction of the lower dimensional operators. 
( See Refs. |3|jl|l for (8L,8fl) operators ). For 
simplicity we set ttik = — M and consider 
the amplitude in the infinite volume. In this case 
the decay amplitude in the method (1) : A'^^^ and 
K ^ TT amplitude in the method (2) : A^^^^ are 
given by 

= M^/f -[a + aC- A(M) + /3(A) • M^] ,(3) 

= M^/f -[a + aD- X{M) + 7(A) • M^] ,(4) 

in one-loop order of CHPT, where / is the pion 
decay constant, A(M) = ■ \og{M^/A^), and 
A is a momentum cut-off for the regularization of 
CHPT. The amplitudes for a finite volume both in 
quenched and full QCD are presented in Refs. (l5[ 

The physical decay amplitude is given by 

A = {ml~ml)/f-[a+aE-\{mK)+S{A)-ml],{5) 

where rriK and ttItt are physical K meson and pion 
mass, is neglected in the bracket in (|^); if this 
is not made, the structure is very complicated, 
which are presented in Refs. [ p^p9| , pOt . 

In (^H) a is a universal constant and inde- 
pendent of kinematics. The term proportional to 
X{M) is the "chiral logarithm term" , and the as- 
sociated constants C, D, and E are determined 
by calculation in one-loop order of CHPT. The 
third terms are the counter terms for the one- 
loop order, which come from tree order diagram 
of the next order of the CHPT operators. The 
associated constants f3, 7, and S depend on the 
kinematics and the momentum cut-off A. 

In usual lattice calculations, the universal con- 
stant a is obtained by a chiral extrapolation of 
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the A^^'' or A^'^'' with the fitting function (||) or 
, where /? and 7 are unknown constants of the 
fitting. Since the precise value of the counter term 
5 in the physical amphtude is not known, it is ne- 
glected in usual lattice calculations. This is the 
cause of a large ambiguity for the lattice predic- 
tion of the amplitude. Determination of counter 
terms from lattice calculations of several ampli- 
tudes are discussed in Refs. |21 2^,^. 

The chiral logarithm term in the amplitude on 
the lattice also depends on the lattice volume and 
whether calculated in quenched or full QCD. Very 
recently, Lin et.al. completed a CHPT calcula- 
tion of the A/ — 3/2 decay amplitudes for Qi, 
Q-j, and Qg operators at general kinematics both 



in quenched and full QCD |23 



4. Recent lattice results 

4.1. Calculation of X ^ 7r(0)7r(0) 

JLQCD collaboration calculated ReA2 from 
the K 7r(0)7r(0) decay amplitude ||. They 
used the plaquette gauge and Wilson fermion ac- 
tion at /3 = 6.1 (1/a = 2.67 GeV) at ttik = Wti- = 
M in the quenched approximation on 24^^ and 32^ 
lattices. They reconstructed the physical decay 
amplitude A^^' from those on the lattice A'"'**- 
by one-loop order of CHPT relation obtained by 
Golterman and Leung p5| It is given by Al^^- ~ 
A^''^-{mK,E) ■G-[1 + M^/{1&tt'^P)-F{EL/2)], 
where rriK is the K meson mass and E is the two- 
pion energy on the lattice, which equals rriK — M 
and E = 2M in their calculations. G is the CHPT 
correction factor in the infinite volume and the 
finite volume correction is given by F{ML) = 
17.827/{ML) + 127rV(AfL)3. 

Very recently, Lin et.al. found a misinterpre- 
tation of the 0(l/i^) term in F{ML) ||. The 
amplitude obtained by the lattice calculation is 
A^ {mK,E + AE), and not A^ {mK,E), where 
AE is the energy shift due to the two-pion in- 
teraction on the finite volume which is given by 
AE = l/{2pL^) in one-loop order of CHPT. 
Thus, in order to use their relation to obtain A^^-, 
we should correct the 0{1/L^) term in F{ML) 
to F{ML) = 17.827/(ML) -I- 2t:^/{ML)^, which 
is same as the correction factor given by Lellouch 
and Liischer (see Sec.^). JLQCD did not consider 
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Figure 1. JLQCD results of Rev42- 

this, but the effect of modification is negligible in 
their simulation points. 

JLQCD results of ReAz • V3 / [G pV^^Vud] ob- 
tained by tree-level (CHPT Tree) and one-loop 
order of CHPT relation (CHPT 1 Loop ) are plot- 
ted in Fig. ^ Here the counter terms of CHPT 
for both the lattice and the physical amplitude 
are neglected. We find that the volume depen- 
dence seen with the tree-level analysis is removed 
after the finite volume corrections at the one-loop 
level. At the same time, the amplitude decreases 
by 30 — 40%. Another noteworthy feature is that 
a sizable lattice meson mass M dependence still 
remains in the amplitude. 

Their final result is ReAa • / [G FV*^Vud] = 
8.9(1.7) X 10-3 GeV^ - 11.4(1.5) x lO'^ GeV^ 
depending on the choice of the scale for the opera- 
tor matching and the CHPT cut-off. These values 
are obtained by a chiral extrapolation of the data 
in Fig. 1^, assuming that the remaining mass de- 
pendence comes only from neglecting the counter 
terms of CHPT on the lattice. JLQCD results are 
consistent with the experiment 10.4 x lO^^^ GeV"^. 
However, a sizable one-loop correction of CHPT 
raises the question whether ignoring higher order 
corrections can be justified. 

4.2. Calculation of iiT — > tt 

The K ^ TT amplitudes have been calculated 
by Pekurovsky and Kilcup ||, RBC [|6|, and 
CP-PACS collaboration [|7). Pekurovsky and 
Kilcup worked with the plaquette gauge and the 
staggered fermions action at f3 — 6.0 (1/a = 
2.1 GeV) and /3 = 6.2 (1/a = 2.8 GeV) in 
quenched QCD, and (3 = 5.7 (1/a = 2.0 GeV) 
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Figure 2. Results for Re^o and Re^2 obtained 
from K ^ TT amplitude. 

for Nf ^ 2 full QCD. RBC and CP-PACS col- 
laboration employed the domain wall fermions in 
the quenched approximation. The gauge action 
is different (RBC : plaquette action, CP-PACS 
: RG-improved action), but the lattice cut-off is 
similar (1/a ~ 2.0 GeV). The three groups set 
rriK = ni-r: = M for the K n amplitude. 

4.2.1. Results for Re^o and Re^2 

The results for Re^o and RcA2 obtained by the 
three groups are summarized in Fig. ^ ; PK stands 
for Pekurovsky and Kilcup, and "Tree" and "1 
loop" refer to the order of CHPT for the phys- 
ical decay amplitude. We should compare the 
results in the same order of CHPT, because it is 
independent of the lattice calculations. The lat- 
tice results for Re^o are almost consistent within 
the three groups, but smaller than the experi- 
ment ReAg = 33.3 x 10~® GeV. It is also seen 
that the results for KeA2 are inconsistent within 
the lattice results. In particular the difference 
between RBC and CP-PACS calculated with the 
same fermion action is about 7a. 

The dominant operators for Re^2 are the Qi 
and • At tuk = TOtt = M, the K ir ampli- 
tudes for these operators are related to the value 
oiBK by (^+|Qf/')|A'+) = (7r+|Q(^/')|X+) = 
Bk ■ 9/4 • [MfKf, where Qf'^^ is the A/ = 3/2 
part of the operator Qj. Thus, ReA2 is propor- 
tional to Bk in the chiral limit approximately. 
In Fig. ^ we show the M dependence for the bare 
Bk obtained from K —> n amplitude by RBC and 
CP-PACS. Here RBC results are those obtained 
by the wall-wall normalization which is used for 
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Figure 3. Bare Bk obtained from K ^ t: ampli- 
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the estimation of KeA2 . 

We attempt to fit the Bk data with the func- 
tion predicted by CHPT : = ail-C-L{M))+ 
/3-M2, where L(M) = l/(167r2/2) . Af^ . log(Af2)^ 

The coefficient of the chiral logarithm term C 
is dealt with in two ways : fixed constant at 
the CHPT prediction (7 = 6 (Fixed logarithm) 
or an unknown constant (Unknown logarithm). 
The CP-PACS data at the largest M, enclosed 
by box in the figure, is omitted in the fitting. 
In Fig. ^ the results of the fitting are tabulated. 
The fitting curves are also plotted by solid lines 
for the "Fixed logarithm" fitting and broken lines 
for the "Unknown logarithm" fitting. The M de- 
pendence of the RBC data are consistent with 
the CHPT prediction, but those of the CP-PACS 
are not. The value of the chiral logarithm given 
by "Unknown logarithm" fitting is inconsistent 
with the CHPT prediction (7 = 6. This inconsis- 
tency with CHPT is also found in previous CP- 
PACS work of Bk In their calculations of 
the K TTTT decay amplitudes, the final results 
are evaluated by the chiral extrapolation with the 
quadratic polynomial function and the "Unknown 
logarithm", which are refereed to as "CP-PACS 
quad." and "CP-PACS log." in Fig. |, respec- 
tively. The final results of RBC are given by the 
"Fixed logarithm" fitting. 
This inconsistency of the M dependence be- 
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Table 1. Results of e'/e ( xlG-* ). 



tween the two groups is reflected in the large dis- 
crepancy of the final results of KeA2. It should 
be noted that the renormalization factor can not 
change the M dependence. A possible reason 
for of this is the lattice cut-off error because the 
gauge actions adopted by the two groups are dif- 
ferent. Numerical investigations on the mass de- 
pendence of Bfc closer to the continuum limit 
are necessary to consolidate calculations of the 
K — > TTTT decay amplitude. 

4.2.2. Results for e'/e 

We found an inconsistency within the lattice 
results for Re^2- The lattice results are also 
troublesome for e'/e. As shown in Table. 1 they 
are far from the experimental values obtained by 
KTeV g and NA48 groups 

One of possible reasons for the discrepancy is 
the quenched approximation. The dominant op- 
erators for e'/e are Qe and Qs- Golterman and 
Pallante pointed out that the naive relation be- 
tween the K —f TT and the K — > tttt decay am- 
plitude for Qe operator shown in (||) and js]) 
are broken in the quenched QCD at 0{M^) [ pl[ . 
As they discussed, this problem can be avoided 
by removing the contraction between q and q 
in Qe = J2q=u,d,Asadb)L^{qbqa)R^- (An alter- 
native method is discussed in Ref. Q). Bhat- 
tacharya et.al. investigated the effect of removing 
the contractions and found the effect to be very 
large for the Qe operator (about 20- 100%) [||. 
CP-PACS collaboration also examined the ef- 
fect for e'/e. They find that while the effect 
is very large, the value of e'/e is still negative 
e'/e = -1.70(53) x 10"^ [||. 

Another possible reason is the lack of final state 
interactions in the K ^ n amplitude, which are 
considered to be important for the A/ =1/2 pro- 
cess. Exploring methods for a direct calculation 
of the decay amplitude from the K — + tttt Green's 
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function is strongly desirable. 

4.3. Calculation at SPQR kinematics 

SPQCDR collaboration calculated the X — > tttt 
decay amplitude from the K tttt Green's func- 
tion at SPQR kinematics directly ||l^ , ^ . In this 
section we show their preliminary results |]35| . 

They calculate the AJ = 3/2 ampfitude for Qi, 
Qt, and Qs which are the dominant operators for 
ReA2 and the A/ = 3/2 part of e'/e. The calcu- 
lations are carried out with the plaquette gauge 
and the Clover- Wilson fermion action with non- 
perturbatively determined Csw at /3 = 6.0 on a 
24^ X 64 lattice. The amplitudes for 7r{p)n{0) at 
— and (27r/i)^ are converted to the physi- 
cal one by the one-loop order of CHPT relations 
obtained in Ref. Finite volume effects have 
not been included in the preliminary results. 

In Fig. ^the ratio between the amplitude M/^ = 
(Tr+Tr^lQilif"*") and the kinematic constant of the 
lowest order of CHPT {{M^+E^)- Mk /2 + A'U ■ 
Ettt ) is plotted. Here M.^ and Mk are pion and K 
meson masses, and E^^ = \/ M"^ The figure 

shows that the ratio depends on the kinematics. 
This means that higher order effects of CHPT 
are significant. Their fit results with the one- loop 
order of CHPT relation are also plotted by open 
symbols. 

The preliminary results for the physical ampli- 
tudes are M4 = 0.0135(80) GeV^ (the experiment 
= 0.0104 GeV^), (Q7(2 GeV)) = 0.14(2) GeV^ 
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and (g7(2 GeV)) = 0.71(6) GeV^ They also 
found that the effects of the one-loop order of 
CHPT are large, 65% for Qi, -34% for Q7, and 
—24% for Qs in their kinematic range. A sizable 
one-loop correction raises the question whether 
ignoring higher order corrections can be justified. 

5. New ideas to solve MT-1 

5.1. Diagonalization method 

The problem of MT-1 also appears in the 
pion 4-point function given by gnm{t) = 
(0|(7r7r)„(t) (7r7r)m(0)|0), where the same defini- 
tion as in (|^) is used. For large time regions 
i 3> the 4-point function behaves as gnm{t) = 
• A, (t) ■ T/^, where = (0|(w)n | (^,) , 
Aj{t) = exp{—Ej ■ t), and |(7r7r)^.) is the energy 
eigenstate with energy Ej. We realize that the 
pion 4-point function contains many exponential 
terms similar to the K — > tttt Green's function. 

This problem can be solved by a diagonal- 
ization of the matrix M(i, to) = 9^^^^{ta) git) 
g~'^^'^{to) at each i, where is some reference 
time and the momenta n and m are regarded as 
matrix indices. The eigenvalues of M{t,to) take 
the form Xj{t,to) = ex.p{—Ej ■ [t — to]), and Vnj = 
(0|(7r7r)„|(^^.) = [g-'/Hto) UHo) A-y^{to)U. 
where U{t()) is an orthogonal transformation ma- 
trix for the diagonalization of M{t, to). 

If the matrix V can be determined precisely 
from the study of the two-pion system, we can 
solve the MT-1 for the — > tttt Green's func- 
tion by defining a new function Gn{t) = X^jKij^ ' 
Gj{t) = An ■ An{t) ■ G^{tK), where Gj{t) is the 
K TTTT Green's function defined by (|^), and 
= ((^^)„|C'|i^). We expect that the new 
Green's function G„(t) behaves as a single ex- 
ponential function for large times and the decay 
amplitude An can be extracted from it easily. 

This diagonalization method was proposed by 
Liischer and Wolf It has been applied to 

many statistical systems ||3^ , |37| , ^ and also to 
the 1 — 2 two-pion system in QCD by Fiebig 
et.al. ^ and by CP-PACS collaboration In 
particular CP-PACS evaluated the 1 = 2 pion 
phase shift with small statistical errors, which is 
extracted from the energy Ej obtained by the lat- 
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AEi = 29.8(22) MeV 

1.000(7) 0.046(6) 
-0.0160(3) 0.324(3) 



L-* = 32^ = (3.2 fm)^ 
p = 380 GeV 

A£i =15.11(50) MeV 

1.000(4) 0.035(2) 
-0.0131(3) 0.352(2) 



= 48' = (4.8 fm)= 
p = 250 GeV 

A£| =5.10(24) MeV 

1.000(6) 0.024(2) 
-0.0094(4) 0.381(3) 




Figure 5. CP-PACS results of two ratio R„{t) 
and Dn{t) for n = 1. 

tice simulation using Liischer's quantization con- 
dition. 

The CP-PACS calculations are carried out with 
the plaquette gauge and the Wilson fermion ac- 
tion in the quenched approximation at /3 = 5.9 
(1/a = 1.9 GeV) on three lattice volumes (24^, 
32^, and 48^). Since the entire matrix of the 4- 
point function gnm{t) cannot be obtained in ac- 
tual lattice calculations, they set a momentum 
cut-off p^jjt = {2-k/LY ■ N and obtain the eigen- 
values \n(t) for n < N. The cut-off dependence 
is also investigated. 

In order to examine the effects of diagonal- 
ization, it is convenient to consider two ra- 
tios, Rn{t) = gnn{t)/[gUt)? and Dnit) = 

An(i, io)/[5n(0]^j where gZ{t) is the pion 2-point 
function with momentum Pn = (27r/L)^ • n. If 
the 4-point function contains only a single expo- 
nential term, then Rn{t) = exp[—AEn ■ t], where 
AEn = En — En and En is twice the n-th pion 
energy. If the momentum cut-off is large enough, 
then Dnit) = exp[~AEn ■ (t ~ to)]- 

In Fig. H the two ratios Rn{t) and Dn{t) for 
n = 1 at m-^jmp = 0.692 are plotted. The pion 
source is located at t = 8. The momentum cut- 
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off is set at TV = 1 and N — 2. It is found that 
the diagonahzation is effective for the small lat- 
tice volume while it is not for the large volumes in 
the figure. The momentum cut-off dependence is 
negligible, however. A single exponential behav- 
ior can be seen for the ratio after diagonahzation 
Dn{t) for all volumes. The values of the energy 
shift Ai5„ and Vnj = (0|(7r7r)„|(7r7r) ■) for the cut- 
off TV = 1 are also tabulated in Fig. That these 
are obtained with small statistical errors demon- 
strate that the MT-1 can be solved by the diago- 
nahzation method. Furthermore their results for 
the 48^ lattice show that the MT-1 is not serious 
for n = 1 on this lattice. Of course the CP-PACS 
investigation is limited to the 1 = 2 two-pion sys- 
tem; the study of the / = system is necessary 
for the A/ =1/2 decay process. 

Finally we comment on the quenched approx- 
imation. It is known from quenched CHPT that 
serious problems appear in the amplitude for the 
two-pion system in quenched QCD due to lack of 
unitarity |^,^ : Chiral and un-physical thresh- 
old divergences appear for both / = and 1 = 2 
two-pion systems, and enhanced finite volume ef- 
fects for 7 = are present. These problems are 
also predicted to appear in the K — s- tttt decay 
amplitudes (l|,|2|,H,|4) . While the presence such 
pathologies has not been numerically confirmed 
in actual lattice simulations, we should make our 
study in full QCD in order to avoid the uncon- 
trollable quenching problems. 

5.2. G-periodic boundary condition 

At this conference Kim ]45| proposed a new 
idea of imposing the G-periodic boundary condi- 
tion for u and d quark and G-periodic boundary 
condition for s quark in the z-direction. The ex- 
plicit forms are given by u{x+z-L) = —i-C(F{x), 
d{x + z ■ L) = +i ■ Cu^{x), and s{x + z ■ L) = 
Cs'^{x), where G is the charge conjugation matrix 
{C~^j^C = 7J). The periodic boundary condi- 
tion is imposed in other directions for all quarks. 
Since pion is an odd eigenstate of G-parity, the 
boundary condition for pion in the z-direction is 
anti-periodic and allowed z-components of mo- 
mentum are restricted to Pz = n/L ■ (2j + 1) 

(J=0,1,---)- 
Kim proposed to extract the decay ampli- 



tude from G{t) = {0\{Tnr){t)O{0)K{tK)\0) , where 
(7r7r)(t) = 'K(p,t)T:{—p,t) and p = ztt/L. Since 
this operator is the interpolating operator for 
the lowest two-pion state under the G-periodic 
boundary condition, G{t) behave as a single ex- 
ponential function for large and the decay am- 
plitude should be extracted easily. Numerical im- 
plementation of this method is future work. 

6. The MT-2 

We finally consider the second part of the 
Maiani-Test problem (the MT-2). As we ex- 
plained above, the MT-1 can be avoided or solved 
by appropriate methods. However, even if we suc- 
ceed in the extraction of the K — > tttt decay am- 
phtude AL'^t- = ((tttt) 1 01 if) from the if -> tttt 
Green's function by lattice simulations, another 
problem remains. The lattice amplitude A^^^' is 
that for the two-pion energy eigenstate on a finite 
Euclidean space-time, and not in the infinite vol- 
ume Minkowski space-time. We can relate A^'^^' 
to the physical amplitude A^^' by using some ef- 
fective theory, but using such theories is the cause 
of large uncertainties in the lattice prediction of 
the decay amplitude as we discussed in the previ- 
ous sections. Recently Lellouch and Liischer ob- 
tained a relation between A^'^^- and A^^- at two- 
pion energy E = tjik Here we show a brief 
derivation. 

The two-pion energy E on the lattice satisfies 
the Liischer quantization condition . Introduc- 
ing a "weak" interaction = nnK ■ A^"^*- , the 
energy is further shifted from E to E' on finite 
volumes, which also satisfies the Liischer quanti- 
zation condition, i.e. 

S{E) = -(j){E) , 5'{E') = -ipiE') mod tt , (6) 

where 5{E) is the pion phase shift for the strong 
interaction, and 5'{E') is that for the strong and 
the weak interactions. The spherical zeta func- 
tion (j) is common for both cases. 

The energy shift E' — E can be estimated by 
perturbation theory for the weak interaction on a 
finite volume. In case oi E = niK-, one finds 

E' -E = ±\A^''^-{E)\ , (7) 

where it should be noted that the energy shift is 
first order in the weak interaction. This is because 
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the energy is degenerate at = 0. We can also 
estimate the difference of the phase shifts at E' by 
perturbation theory in the infinite volume, which 
yields 

6'{E') - 5{E') = — . ^ . ' , , (8) 

^ ' ^ ' Sir E' P^ + mj^ ' ^ ' 

where p' = \J E'"^ /A — m'^, and the denomina- 
tor is the S'-channel K meson propagator with 4- 
momentum Pr- = {iE',0). Substituting (||) into 
the Liischer quantization condition and ex- 
panding for E' — E hy using (|^), we obtain the 
following formula at iJ = rriK (LL-formula) . 

|APh-(^)|2 = c{E) ■ p{E) ■ , (9) 

where C{mK) — ?>2t:'^ ■ m\/p, p = \/ E'^ /4, — mj, 
and p{E) = E/{AttP) ■ A/{Ap)[5{E) + (j){E)]. 

The validity of the LL-formula for the A/ = 
3/2 amplitude is confirmed in one- loop order of 
CHPT even in the quenched theory ||2^. This is 
not the case for the A/ —1/2 amplitude. Since 
the correction factor depends on the weak opera- 
tors due to the lack of unitarity in quenched QCD, 
the simple relation from A^'^^- to A^^- can not be 
obtained HQQ. 

It should be noted that the LL-formula given by 
is relied on the Liischer quantization condition 
which holds for the energy on a finite periodic 
box in the two-pion center of mass frame. The 
formula is not applicable to the SPQR kinematics 
which is in the laboratory frame. An extension 
of the LL-formula to the general frame can be 
obtained easily using the Rummukaine- Gottlieb 
quantization condition for the general frame . 

We note that the formula can not be applied 
for the G-periodic boundary condition, because 
the anti-periodic boundary is imposed in the z- 
direction. Extending the formula to this case is 
not trivial, since the condition on the energy in 
such a boundary condition is not known. 

Lin et al. derived the LL-formula from a dif- 
ferent approach Q. They found the following 
simple relation between the two-pion state \-k-k) 
in the physical decay amplitude and that on the 
lattice Ktttt)) at general energy E : Ktttt)) 
(47r) • ^p{E) ■ E/p ■ |(^). This gives the exten- 
sion of the LL-formula to the case oi E ^ ttik- 



They also investigated the validity of the LL- 
formula, and found that the volume has to be 
sufficiently large or the two-pion interaction has 
to be weak enough, so that the boundary condi- 
tion does not distort the two-pion wave function. 
This is also required for the Liischer quantization 
condition. 

At present there are few studies of the two-pion 
system, and informations are scarce of the two- 
pion wave functions. It is extremely important to 
investigate the two-pion system, before embark- 
ing on a numerical application of the LL-formula 
to the K TTTT decay amplitude. 

7. Summary 

In this article we discussed recent theoretical 
and numerical progresses in the calculation of the 
K ^ TTTT decay amplitude. While this calculation 
has been plagued with a number of difficulties, we 
discussed that most of them had been shown to 
be theoretically solvable. In particular, the MT-1 
can be avoided by a judicious choice of kinemat- 
ics and boundary conditions or can be solved by 
diagonalizations, in principle. 

In our opinion, the most desirable way of cal- 
culating the decay amplitude is the following. 
We extract the decay amplitude from the K —> 
TTTT Green's function directly using the SPQCDR 
method or one of new ideas described in previous 
sections. We then reconstruct the physical decay 
amplitude from that on the lattice by the LL- 
formula or its extensions without using effective 
theories of QCD. 

In this procedure, a quantitative study of the 
two-pion system is necessary to examine the va- 
lidity of the LL-formula. Also all calculations 
should be made in full QCD to avoid uncontrol- 
lable quenching problems due to the lack of uni- 
tarity in the quenched QCD. 

The final obstacle to avoid using effective the- 
ories in the actual lattice calculations is chiral 
extrapolation. Practical lattice simulations are 
carried at - (400 - 800) MeV. We have 
to extrapolate our decay amplitude obtained by 
the LL-formula to that at the physical mass with 
some assumption for the mass dependence. Here 
effective theories such as CHPT is used. This 
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causes large uncertainties in the lattice predic- 
tion of the decay amplitude, because there is no 
reliable effective theory for such a long extrap- 
olation. A unique solution to this problem is a 
simulation at or near the physical point. 
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